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The  93rd  Meeting  of  the  Acoustical  Society  of  America  Invited  Paper  KK3, 

2 P.M. , June  9,  1977.  Directional  Patterns  of  Transducer  Arrays. 

Paul  M.  Kendig,  Applied  Research  Laboratory,  Pennsylvania  State  University, 
University  Park,  Penna.  16802 

In  order  to  locate  objects  under  water  by  means  of  sound  waves,  it  is 
desirable  to  transmit  sound  in  - or  receive  sound  from  - one  principal 
direction.  This  is  usually  accomplished  by  employing  an  array  of  electro- 
acoustic transducer  elements.  If  the  object  to  be  located  is  itself  a 
source  of  sound,  then  only  a receiving  array  (hydrophone)  is  required.  If 
the  object  is  not  a self  radiator  of  sound,  then  the  array  must  transmit  a 
sound  pulse  and  determine  the  direction  of  the  object  from  the  direction 
of  the  returned  echo.  For  this  (the  active  case)  the  same  array  usually 
serves  both  functions,  i.e.,  sound  projector  and  receiving  hydrophone.  The 
range  of  the  object  may  also  be  determined  with  the  active  system  by 
measuring  the  travel  time  of  the  sound  pulse  to  the  target  and  back. 

One  characteristic  of  a reciprocal  electroacoustic  transducer  (the 
type  most  commonly  used)  is  that  the  directional  responses  or  patterns  are 
the  same  for  both  the  transmitting  and  receiving  modes.  This  is  convenient 
because  in  the  following  discussion,  it  is  not  necessary  to  distinguish 
between  the  two  modes  because  every  thing  that  follows  applies  to  both. 

The  following  discussion  is  given  in  greater  detail  in  Reference  1. 

A typical  arrangement  of  transducer  elements  in  a planar  array  is 
shown  in  Figure  1 where  the  small  circles  represent  the  locations  of  the 
centers  of  the  radiating  faces  of  the  transducing  elements.  These  are 
usually  square  or  rectangular  surfaces.  The  spacing  of  columns  and  rows 
are  uniformly  spaced  but  not  necessarily  equal  to  each  other. 

The  directional  pattern  in  a plane  containing  the  :<-axis  is  Identical 
to  that  of  a line  parallel  to  the  >:-axis  and  with  strengths  equal  to  the 
sums  of  those  in  the  respective  columns.  We  can  tlicn  determine  this 
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specific  pattern  by  reducing  the  problem  to  that  of  determining  the  pattern 
of  a line  of  equally  spaced  sources.  Of  course,  the  same  applies  to  a 
pattern  in  a plane  that  contains  the  y-axls. 

A slight  modification  to  the  pattern  obtained  by  treating  the  elements 
as  point  sources  is  due  to  the  directional  characteristics  of  the  radiating 
faces  of  the  elements  themselves  which  are  all  assumed  to  be  identical. 

The  actual  pattern  is  then  the  product  of  the  pattern  of  one  of  the  identical 
elements  and  the  pattern  of  the  line  array  of  point  sources  (first  product 
theorem) . 

For  the  same  reasons  that  the  pattern  of  the  planar  array  in  a plane 
parallel  to  either  lines  or  rows  reduce  to  that  of  a row  of  sources,  the 
pattern  of  a rectangular  source  in  a plane  parallel  to  a side  reduces  to 
that  of  a line.  Consequently,  we  shall  first  determine  the  pattern  of  a 
line  source.  Parenthetically,  one  should  be  reminded  that  directional 
patterns  are  simply  diffraction  patterns  and  as  such  are  similar  to  much 
that  we  have  learned  in  optics. 

Figure  2 helps  describe  the  terms  used  to  determine  the  pattern  of  a 
line  source  in  a plane  that  contains  the  line.  It  shows  a plane  wave 
approaching  the  line  at  an  angle  of  0 with  the  normal  to  the  line.  The 
signal  received  at  the  incremental  portion  of  the  line  dx  leads  the 
signal  (reference  phase)  received  at  the  origin  by  an  angle  4/  which  is 
shown  to  be 

sin  0 = kx  sin  0 

where  X is  the  wavelength  of  the  sound  wave  and  x is  the  distance  of 
the  element  measured  from  the  origin. 
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The  pattern  response  is  given  by 


P^  = Q 


i k X sin  9 


where  Q is  the  source  strength  per  unit  length  and  i is  the  length  of 
the  line.  Integrating,  we  get 


"a  = ^ 


i kH/2  sin  9 -i  kZ/2  sin  9 


i k sin  9 


which,  from  an  identity  expressing  the  sine  of  an  angle  in  terms  of  exponentials. 


becomes 


. ki-  . Q 
sin  ~ sin  9 


k£  . ^ „ Ql  sin  u'  ^ 

Let  u = sin  0 . Then  P^  = (17 

which  is  just  the  same  as  the  diffraction  of  a slit  in  optics. 

An  important  relationship  for  a plane  rectangular  radiating  surface  is 

obtained  in  the  case  where  the  source  strength  F(x,y)  = f(x)  g(y)  . In 

this  case 


i k(x  cos  a + y cos  B) 


P^  = F(x,y)  e 


dx  dy 


i k X cos  a -ire  iky  cos  6 

) e dx  g(y)  e dy  (2) 


where  a , B , v are  angles  between  direction  of  observation  and  the 


X,  y and  z <txos,  rcsju'ct  ive! y. 
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This  very  important  theorem  states  that  if  the  source  strength  for  a 
plane  rectangular  surface  may  be  expressed  as  the  product  of  a function 

f(x)  and  a function  f(y)  , then  the  directional  pattern  is  the  product 

of  two  patterns  in  planes  containing  the  normal  to  the  surface  and  parallel 
to  the  X and  y axes,  respectively.  The  full  significance  of  this  theorem 
will  appear  later  when  we  discuss  rectangular  arrays  of  elements. 

Since  we  have  already  shown  that  patterns  in  the  principal  planes  of 
a rectangular  array  of  point  sources  can  each  be  given  by  a line  of  point 
sources,  we  shall  consider  an  equally  spaced,  in-phase  line  of  sources. 

The  symbols  that  shall  be  used  for  such  arrays  are  shown  in  Figure  3.  The  i 

incident  sound  wave  meets  the  line  array  at  the  angle  0 as  shown.  All  I 

phases  shall  be  referenced  to  that  of  the  signal  received  by  an  element  at  j 

the  origin  whether  one  is  actually  there  or  not.  Note  the  right-left  | 

symmetry  of  the  signal  amplitudes. 

Figure  4 is  a representation  of  the  signals  recleved  by  each  of  the 
elements  in  a six  element  line  array  where  u = sin  6 . The  vectors 

above  the  axis  represent  the  signals  received  from  the  elements  on  the 
right  and  those  below  the  axis,  the  signals  from  the  left. 

It  is  easily  seen  from  Figure  4 that  the  sum  pattern  of  the  six  element 
line  source  is  given  by,  i 

P , = A,  cos  u + A„  cos  3u  + A_  cos  5u  (3) 

00  1 2 3 

and  that  over  the  major  lobe  the  phase  of  the  entire  line  array  is  the 
same  as^that  of  the  reference  signal.  As  0 and  consequently  u Increases, 

P^^  will  at  one  point  become  ?,ero  (a  null  in  the  pattern)  and  then  the 
phase  suddenly  reverses  to  180°.  Further  increases  of  0 and  u produce 


more  nulls  and  phase  reversals. 
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Figure  4 also  shows  that  the  difference  pattern  is  given  by, 

Pp  , = A,  sin  u + A„  sin  3u  + A.  sin  5u 

0,0  1 Z J 

In  this  case,  we  note  that  for  small  angles  of  incidence  on  the  right  the 
phase  leads  by  90°  and  on  the  left  lags  by  90°.  Whereas  the  sum  pattern 
is  a maximum  for  0=0,  the  difference  pattern  is  a null. 

Corresponding  pattern  functions  for  odd  number  of  elements  are: 

P = A„  + A,  cos  2u  + A-  cos  4u  + . . . 

On  0 1 2 


P.  = A,  sin  2u  + A-  sin  4u  + . . . 
on  1 2 

The  characteristics  of  P , and  P-  . are  illustrated  in  Figure  5.  The 
combined  use  of  these  patterns  give  the  necessary  information  for  locating 
underwater  targets.  Note  that  over  a range  of  incident,  angles  near  0°, 
the  difference  between  the  magnitudes  of  the  sum  and  difference  patterns 
increases  nearly  linearly  with  0 . If  the  difference  pattern  leads  the 
sum  pattern  by  90°,  the  target  is  on  the  right  and  if  it  lags  by  90°,  the 
target  is  on  the  left.  In  practice,  the  difference  pattern  is  usually  phase 
shifted  90°  so  that  it  is  in  phase  for  targets  on  one  side  and  out  of  phase 
for  the  other. 

It  must  be  emphasized  that  those  phase  relationships  hold  only  over  a 
relatively  small  range  of  incident  angles  on  both  sides  of  the  normal.  In 
the  min(^  lobe  regions  we  noted  earlier  that  phase  reversals  occur  as  0 is 
varied  and  where  they  occur  are  not  the  same  for  both  patterns.  Consequently, 
minor  lobes,  unless  drastically  reduced,  will  give  false  information.  For 
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these  reasons  and  others,  methods  for  reducing  minor  lobes  have  been  developed 
which  are  known  as  tapering  or  shading  the  array.  It  consists  of  a systematic 
variation  of  the  signal  amplitude  coefficients  of  the  elements.  The  process 
will  be  illustrated  by  an  application  to  the  six  element  array  illustrated 
above . 

Note  that  the  pattern  function  was  given  by 


P , = A,  cos  u + A-  cos  3u  + A_  cos  5u 
06  1 2 3 


(3) 


which  is  a Fourier  series.  This  suggested  the  use  of  a Fourier  transform 
relating  the  source  strengths  of  the  array  to  the  pattern 
function.  In  one  application,  an  amplitude  distribution  was  derived  to 
provide  the  Gaussian  pattern.  In  this  case,  the  Fourier  transform  is  Itself 
a Gaussian  function. 

However,  the  most  common  method  of  tapering  is  that  developed  by 
Dolph  (Reference  2)  and  generally  known  as  the  Dolph-Cbebyshev  method  because 
it  sets  the  pattern  function  above  equal  to  an  appropriate  Chebyshev  poly- 
nomial . 

In  Equation  (3)  above,  we  let  x = cos  u and  then  express  the  cosine 
functions  in  powers  of  x using  Chebyshev  polynomials.  Thus,  Equation  (3) 
becomes 


P , = A,x  + A (4x^  - 3x)  + A (16x^  - 20x^  + 5x) 

0,6  12  3 

= 16A^x^  + (4A2  - 20A2>x^  + (A^  - 3A^  + 5A^)x  . (4) 

The  use  of  Chebyshev  pol  ynom  I .i  1 s abo’.'e  Is  somewhat  trivial  comj'nrcd  to  the 
really  significant  use  whicli  is  to  set  the  array  polynomial  above  equal  to 
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the  Chebyshev  polynomial 


T, (x  x)  = 16x  ^x^  - 20x  ^x^  + 5x  X 

5 o o o o 


(5) 


where  x >1  and  is  given  by  x = cosh{—  cosh  ^r} 


(6) 


where  r = lobe  ratio, 

minor 


The  amplitudes  and  A^  are  now  easily  calculated  by  equating 

coefficients  of  like  powers  of  Equations  (4)  and  (5). 

Figure  6 Illustrates  the  properties  of  the  fifth  degree  Chebyshev 
polynomial  which  applies  to  the  six  element  array.  Over  the  range  of 
-1  ^ x^x  1 » *^he  Chebyshev  polynomial  function  ranges  between  +1  because 
in  this  range  the  function  is  limited  to  the  cosine  function.  However,  by 
extending  the  independent  variable  x^x  beyond  +1,  the  Chebyshev  polynomials 
can  take  values  to  + infinity.  Thus,  we  can  make  x^x  and  consequently 
T^(x^x)  as  large  as  possible.  T^ (x^x)  will  then  correspond  to  the  peak  or 
maximum  of  our  response.  On  the  other  hand,  the  minor  lobes  will  all  be 
the  same  and  equal  to  unity. 

This  design  method  is  optimum  in  the  sense  that,  for  a given  minor  lobe 
reduction,  one  obtains  the  narrowest  possible  beamwidths  for  the  equally 
spaced,  in-phase  elements  at  a given  frequency.  We  design  so  that  at 
0 = 90°,  the  Chebyshev  polynomial  is  of  unity  magnitude  just  beyond  the 
last  null. 

There  are  relatively  narrow  limits  on  the  range  of  element  spacing. 


If  the  spacing  is  too  large,  e.g.,  greater  than  X , a second  major  lobe 
would  appear  before  0 = 90°,  Also,  the  method  does  not  work  if  the  spacing 
Is  much  less  than  X/2  . For  smaller  spacings,  a slightly  different  meth.od 
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was  introduced  by  Riblet  (Reference  3)  and  extended  by  Pritchard  (Reference  4). 

In  this  method  which  can  be  applied  to  an  odd  number  of  elements,  a Chebyshev 
polynomial  in  terms  of  an  angle  <j>  = sin  6 is  used  Instead  of  u = ^ sin  6 . 

This  analysis  leads  to  narrower  patterns  than  those  given  by  the  Dolph-Chebyshev 
technique  and  are  called  super-directive  arrays.  However,  as  the  spacing 
becomes  smaller,  some  of  the  excitation  amplitudes  become  negative,  that  is, 
out  of  phase  with  the  others.  For  very  small  spacing,  the  sum  is  considerably 
less  than  the  sum  of  absolute  values.  Smaller  tolerances  on  the  amplitudes 
are  then  generally  required  and  over-all  efficiency  is  decreased. 

Some  characteristics  of  the  Dolph-Chebyshev  patterns  can  be  determined 


directly  from  the  values  of  x^  , d and  X . Of  course,  the  latter  depends 
on  the  frequency.  These  characteristics  are  (1)  maximum  allowable  spacing, 
and  (2)  beam  width  at  any  specified  levels  below  the  peak.  These  results 
are  shown  in  Figure  7.  Our  experience  with  transducer  elements  cemented  to 
a rubber  pad  and  mounted  in  a housing  showed  that  since  there  is  a baffle 
effect  that  significantly  reduces  levels  at  incident  angles  beyond  70°,  it 
was  possible  to  design  for  the  last  minor  lobe  to  occur  around  70°.  The 
dashed  line  Indicates  the  extension  of  this  tolerance. 

We  shall  now  Illustrate  the  calculations  of  the  excitation  amplitudes 


for  our  six  element 

array. 

If 

we  design 

for  a 

30  dB  minor 

lobe 

reduction , 

20  log  r will  be  set 

equal 

to 

30,  which 

gives 

r = 31.6  . 

From 

Equation  (6) 

we  get  a value  of  1. 

35  for 

X 

0 

. From  Equations  (4)  and  ( 

5)  we 

obtain  tlie 

following  values  for 

\ • 

= 15,60 

= 

10.70 

^^3  = 

4.66 

normalised  Aj  = 1.000 

"2 

= 

0.685 

A.  = 

0.  298 
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Figure  8 illustrates  how  these  values  are  used  to  extend  the  line  array 
design  to  that  for  a two  dimensional  array  using  Equation  (2).  Due  to  the 
symmetry,  only  one  quadrant  of  the  array  is  necessary.  Figure  9 gives  the 
normalized  values  for  these  elements.  Frequently  the  outside  corner  element 
is  deleted  so  as  to  confirm  to  a circular  housing.  The  pattern  deterioration 
is  generally  small  but  in  any  case  its  effect  is  easily  calculated.  The 
numbers  in  parentheses  give  one  possible  set  of  turns  on  magnetostrictive 
stacks  or  on  the  transformer  secondaries  used  with  electric  coupled  transducer 
elements. 

The  coefficients  that  provide  an  optimum  sum  pattern  are  unsuitable  for 
an  optimum  difference  pattern.  In  fact,  they  degrade  the  difference  pattern. 
Consequently,  the  difference  pattern  must  be  Independently  optimized. 

With  a slight  modification,  the  Dolph  technique  may  be  employed  to 
optimize  a difference  pattern,  i.e.,  where  the  two  halves  are  in  phase 
opposition.  The  difference  pattern  function. 


Pr  , = A,  sin  u + A„  sin  3u  + A.  sin  5u 

0,6  12  3 


may  be  expanded  in  terms  of  sin  u as  follows: 


P.  , A,  sin  u + A„(3  sin  u - A sin^  u)  + A (5  sin  u - 20  sln^  u + 16  sin^ 
w y D X 3 


Now,  if  we  divide  by  sin  u , we  get. 


= A,  + A„(3  - 4 sin^  u)  + A_(5  - 20  sin^  u + 16  sin^  u) 
sin  u 1 z i 


2 2 

If  we  substitute  1 - cos  u for  sin  u and  x tor  cos  u , the  expression 


-10- 
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above  becomes , 

-7-^  = (A,  - A,  + AJ  + (4A„  - 12A.)  + 16A_x^  . 

sin  u 1 2 3 2 3 3 

We  now  equate  this  relationship  to  a 4th  degree  Chebyshev  polynomial  and 

proceed  as  before  with  the  sum  pattern.  Please  recall,  however,  that  we 

used  a 5th  degree  polynomial  for  the  sum  pattern.  Also,  note  that  since 

Pp  , = sin  u T, (x  x)  , it  is  zero  at  0=0  and  then  the  absolute  value 
0,6  4 o 

Increases  rapidly  for  moderately  small  values  of  0 on  both  sides  of  zero. 

The  use  of  Equation  (2)  to  extend  the  line  array  to  a plane  is  shown 
in  Figure  10.  Here,  the  B factors  are  just  the  excitation  coefficients 
for  an  optimum  sum  pattern.  In  fact,  if  all  elements  were  connected  in 
phase,  we  would  obtain  an  optimum  vertical  sum  pattern,  but  when  the  right 
and  left  halves  are  connected  in  phase  opposition,  we  obtain  the  optimum 
horizontal  difference  pattern. 

Figure  11  gives  relative  source  amplitudes  for  the  optimum  difference 
pattern.  The  distinct  difference  between  these  coefficients  and  those  for 
the  sum  pattern  is  quite  evident  by  noting  that  the  elements  in  the  second 
column  now  have  the  greatest  excitation  instead  of  the  first  column  as  in 
the  case  of  the  sum  pattern. 

No  direct  method  of  finding  x was  available  for  the  difference  pattern 

o 

but  minor  lobe  reductions  could  be  calculated  for  assumed  values  of  x 

o 

(see  Reference  1)  and  plotted  as  shown  in  Figure  12.  Thus,  these  plots 
permit  us  to  select  the  appropriate  x^  for  a given  minor  lobe  reduction. 

Figure  5 presents  a linear  plot  of  the  optimum  sum  and  difference 
patterns.  When  this  method  for  the  optimum  difference  pattern  was  developed, 
it  appeared  that  it  would  apply  only  to  an  even  number  of  array  columns. 
However,  Geoffrey  Wilson  (Reference  5)  has  shown  that  it  can  be  applied  to 
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an  odd  number  of  columns,  in  which  case  the  central  column  simply  is  not 
used  to  form  the  difference  pattern. 

Figure  13  is  of  academic  interest.  It  represents  the  locus  of  the 
resultant  signal  from  1/2  of  a transducer  designed  for  (a)  an  optimum  sum 
pattern  and  (b)  an  optimum  difference  pattern.  It  is  the  resultant  of  the 
vectors  shown  in  Figure  4.  The  sum  pattern  response  is  the  projection  on 
the  x-axls  and  that  of  the  difference  pattern  is  the  projection  on  the  y-axis. 
It  clearly  shows  why  the  minor  lobes  are  so  small. 

Tilted  beams  and  end-fire  arrays  are  obtained  by  introducing  delay  lines 
or  phase  shift  networks  into  the  circuits  of  each  element  or  maybe  column 
or  row  of  the  array  so  that  for  a given  specified  direction,  the  outputs 
of  signals  (including  delays  or  phase  shifts)  will  be  in  phase.  The 
tapering  procedures  are  generally  similar  to  those  described  earlier. 

Most  of  the  development  for  end-fire  arrays  (beam  tilted  to  90°) 
resulted  in  superdirective  arrays  with  some  characteristics  similar  to 
those  of  Riblet  (Reference  3)  and  Pritchard  (Reference  A). 

A somewhat  different  approach  to  the  production  of  tilted  sum  and 
difference  patterns  was  developed  by  W.  J.  Hughes  and  W.  Thompson,  Jr. 
(Reference  6).  I shall  demonstrate  its  use  by  an  application  to  a tilted 
sum  pattern  which  may  be  written, 

P^  = cos(u  - (J))  + A2  cos  3(u  -()))+  A^  cos  5(u  - 4>)  + . . . 

where  (p  is  the  required  phase  delay  for  elements  on  one  side  of  center 
and  phase  advance  for  elements  on  the  other  side.  Using  the  trigonometric 
Identity  for  the  cosine  of  the  difference  of  two  angles,  the  equation  above 


becomes , 
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A 


p 

a 


= cos  (fi  cos  u + cos  3({)  cos  3u  + cos  54>  cos  5u  + . . . ] 

+ [A  sin  <J)  sin  u + A^  sin  3<J)  sin  3u  + A^  sin  5ij)  sin  5u  + . . . ] 


Now  for  a given  tilt  angle  (j)  , the  sines  and  cosines  involving  (])  are 
constants.  Therefore,  the  pattern  function  becomes 


COS 

u + 

"2 

COS 

3u 

+ 

"3 

COS 

5u 

+ . 

. .] 

[A^ 

sin 

u + 

^2 

sin 

3u 

+ 

"3 

sin 

5u 

+ . 

. .] 

where  A'  = A cos  6 , and  A"  = A sin  > etc. 
n n n n 

We  note  that  P^  is  now  the  sum  of  a sum  pattern  and  a difference 

pattern.  Since  the  latter  is  in  phase  quadrature  with  the  former,  it  must 

be  shifted  90°  before  combining  the  two.  Although  the  method  is  applicable 

to  any  values  of  A , tapering  is  easily  achieved.  The  unprimed  A 
m n 

coefficients  are  determined  as  before  by  using  Chebyshev  polynomials.  The 
Hughes  - Thompson  technique  is  also  applicable  to  tilted,  difference  patterns 
and  end  fire  arrays.  It  is  interesting  to  note  that  all  pattern  tapering 
techniques  discussed  so  far  have  used  Chebyshev  polynomials. 

An  array  tapering  for  an  entirely  different  purpose  was  developed  by 
W.  J.  Trott  (Reference  7).  Its  purpose  was  to  provide  plane  waves  in  a 
region  of  the  near  field  of  a relatively  large  array  in  order  to  calibrate 
transducers  in  the  near  field. 

It  is  well  known  that  the  intensity  along  the  axis  of  a circular 
radiating  piston  whose  diameter  is  large  compared  to  the  wavelength  has 
numcroi^  maxima  and  minima  (actually  nulls)  in  the  near  field.  It  is  also 
known  that  the  total  number  of  side  lobes  appearing  in  the  far  field  is 
just  twice  the  number  of  maxima  which  occur  directly  in  front  of  the  piston. 
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The  above  remarks  suggest  that  if  minor  lobes  are  eliminated  or  greatly  reduced, 
this  axial  variation  should  likewise  be  eliminated  or  reduced.  Indeed  this 
is  the  case.  Applying  these  principles,  Trott  adjusted  the  excitation 
amplitudes  of  a large  array  of  elements  in  order  to  provide  an  essentially 
plane  wave  region  directly  in  front  of  the  array.  In  his  example,  a 100  cm 
diameter  array  provided  a cylindrical  space  in  front  of  the  array  about  50  cm 
in  diameter  and  about  the  same  in  length  wherein  there  were  plane  waves. 

Although  it  was  not  expressly  stated,  all  the  previous  discussion  was 
based  on  linear  acoustic  waves.  When  very  intense,  so-called  finite  sound 
waves  are  projected,  non-linear  or  distorted  waves  are  produced  with  a whole 
host  of  interesting  properties.  For  example,  if  two  pure  tones  are  projected 
into  the  same  medium,  harmonics  of  the  waves  are  produced,  as  well  as  sum 
and  difference  frequency  waves.  The  highly  directive  character  of  these 
difference  frequency  beams  was  first  recognized  by  Peter  J.  Westervelt 
(Reference  8 ) . 

In  the  following  discussion  we  are  primarily  Interested  in  the  difference 
frequency  that  develops  in  the  projected  acoustic  field.  Usually  the 
frequencies  of  the  high  intensity  primary  waves  are  close  together  so  that 
the  difference  frequency  is  much  smaller,  say  one-tenth  or  less. 

When  the  near  field  primary  wave  absorption  is  quite  small,  the  inter- 
action zone  will  extend  beyond  the  mean  Rayleigh  distance  r^  . The  inter- 
action zone  is  that  region  in  front  of  the  transducer  face  where  wave 
distortion  occurs  and  acoustical  energy  is  pumped  from  the  primary  weaves 
into  the  distortion  waves,  including  the  difference  frequency  wave.  Since 
this  lifter  is  of  much  lower  frequency  than  any  of  the  other  waves  and 
consequently  attenuated  miuh  less  than  the  others,  it  usually  will  be 
transmitted  to  greater  ranges  tlian  the  otliers  even  thougii  its  equividcnc 
field  pressure  is  much  below  that  of  the  jirlmary  waves.  The  so-called 
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pumping  occurs  until  the  spreading  loss  of  the  primary  waves  reduce  the 
Intensity  to  such  a level  that  distortion  no  longer  exists.  Since  the 
virtual  sources  that  produce  the  difference  frequencies  are  essentially 
limited  to  the  main  beams  of  the  primary  waves,  the  low  frequency  difference 
pattern  will  be  essentially  the  same  as  those  of  the  much  higher  frequency  primary 
patterns  or  even  a little  sharper.  In  fact,  its  far  field  directivity 
response  is  essentially  equal  to  the  product  of  the  directivity  responses 
of  the  primary  patterns.  One  very  Important  consequence  of  this  is  the 
almost  total  absence  of  minor  lobes  in  the  difference  frequency  pattern. 

In  order  to  achieve  the  same  directivity  with  a conventional  source 
operating  at  the  difference  frequency,  the  source  diameter  would  have  to 
exceed  the  diameter  of  the  parametric  array  by  the  ratio  of  the  mean 
primary  frequency  to  the  difference  frequency. 

The  bandwidth  of  the  difference  frequencies  may  be  very  broad  because 
it  is  really  determined  by  the  Q of  the  transducer  that  projects  the 
primary  waves. 

Somewhat  similarly,  difference  frequency  waves  can  be  produced  which 
have  a constant  beam  width  over  very  wide  frequency  bands  (e.g.,  greater 
than  two  octaves).  Of  course,  sum,  difference  and  tilted  beams  can  be 
produced . 

When  the  absorption  coefficient  of  the  primary  waves  is  large,  the 
attenuation  may  be  so  great  that  all  of  the  interaction  may  take  place 
within  the  near-field  of  the  transducer,  i.e.,  at  ranges  less  than  r 

o 

In  this  case,  the  interaction  is  essentially  confined  to  a cylinder  whose 
diameter  is  that  of  the  transducer  radiating  face  and  of  length  roughly 
determined  by  the  absorption  coefficient.  This  interaction  zone  is  really 
an  extension  of  the  source  itself.  Indeed,  the  beam-width  of  this  virtual 
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end-fire  array  is  given  by  the  same  relationship  as  that  for  an  end-fire 
array  and  is  always  less  than  that  for  the  "spreadlng-loss-llmited"  case 
Parametric  receiving  arrays  are  formed  in  a fluid  by  projecting  a 
finite-amplitude  "pump  wave"  into  the  medium  to  serve  as  a "carrier"  for 
the  weak  incoming  signal  whose  frequency  equals  the  difference  frequency 
of  the  transmitted  waves. 
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Figure 


Figure  7 - Beam  width  vs.  spacing  for  30  db  minor  lobe  reduction 
Parameter;  number  of  elements. 
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The  93rd  Meeting  of  the  Acoustical  Society  of  America  Invited  Paper  KK3, 

2 P.M. , June  9,  1977.  Directional  Patterns  of  Transducer  Arrays. 

Paul  M.  Kendlg,  Applied  Research  Laboratory,  Pennsylvania  State  University, 
University  Park,  Penna.  16802 

In  order  to  locate  objects  under  water  by  means  of  sound  waves.  It  Is 
desirable  to  transmit  sound  In  - or  receive  sound  from  - one  principal 
direction.  This  is  usually  accomplished  by  employing  an  array  of  electro- 
acoustic transducer  elements.  If  the  object  to  be  located  is  itself  a 
source  of  sound,  then  only  a receiving  array  (hydrophone)  is  required.  If 
the  object  Is  not  a self  radiator  of  sound,  then  the  array  must  transmit  a 
sound  pulse  and  determine  the  direction  of  the  object  from  the  direction 
of  the  returned  echo.  For  this  (the  active  case)  the  same  array  usually 
serves  both  functions,  l.e.,  sound  projector  and  receiving  hydrophone.  The 
range  of  the  object  may  also  be  determined  with  the  active  system  by 
measuring  the  travel  time  of  the  sound  pulse  to  the  target  and  back. 

One  characteristic  of  a reciprocal  electroacoustic  transducer  (the 
type  most  commonly  used)  is  that  the  directional  responses  or  patterns  are 
the  same  for  both  the  transmitting  and  receiving  modes.  This  is  convenient 
because  in  the  following  discussion,  it  is  not  necessary  to  distinguish 
between  the  two  modes  because  every  thing  that  follows  applies  to  both. 

The  following  discussion  is  given  in  greater  detail  in  Reference  1. 

A typical  arrangement  of  transducer  elements  in  a planar  array  is 
shown  in  Figure  1 where  the  small  circles  represent  the  locations  of  the 
centers  of  the  radiating  faces  of  the  transducing  elements.  These  are 
usually  square  or  rectangular  surfaces.  The  spacing  of  columns  and  rows 
are  uniformly  spaced  but  not  necessarily  equal  to  each  other. 

The  directional  pattern  In  a plane  containing  the  ;<-axis  is  identical 
to  that^of  a line  parallel  to  the  >:-axis  and  with  strengths  equal  to  the 
sums  of  those  In  the  respective  columns.  We  can  then  determine  tin's 
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specific  pattern  by  reducing  the  problem  to  that  of  determining  the  pattern 
of  a line  of  equally  spaced  sources.  Of  course,  the  same  applies  to  a 
pattern  in  a plane  that  contains  the  y-axis. 

A slight  modification  to  the  pattern  obtained  by  treating  the  elements 
as  point  sources  is  due  to  the  directional  characteristics  of  the  radiating 
faces  of  the  elements  themselves  which  are  all  assumed  to  be  identical. 

The  actual  pattern  is  then  the  product  of  the  pattern  of  one  of  the  identical 
elements  and  the  pattern  of  the  line  array  of  point  sources  (first  product 
theorem) . 

For  the  same  reasons  that  the  pattern  of  the  planar  array  in  a plane 
parallel  to  either  lines  or  rows  reduce  to  that  of  a row  of  sources,  the 
pattern  of  a rectangular  source  in  a plane  parallel  to  a side  reduces  to 
that  of  a line.  Consequently,  we  shall  first  determine  the  pattern  of  a 
line  source.  Parenthetically,  one  should  be  reminded  that  directional 
patterns  are  simply  diffraction  patterns  and  as  such  are  similar  to  much 
that  we  have  learned  in  optics. 

Figure  2 helps  describe  the  terms  used  to  determine  the  pattern  of  a 
line  source  in  a plane  that  contains  the  line.  It  shows  a plane  wave 
approaching  the  line  at  an  angle  of  0 with  the  normal  to  the  line.  The 
signal  received  at  the  incremental  portion  of  the  line  dx  leads  the 
signal  (reference  phase)  received  at  the  origin  by  an  angle  ip  which  is 
shown  to  be 


ip  = sin  0 = kx  sin  0 

where  A is  the  wavelength  of  the  sound  wave  and  x is  the  distance  of 
the  element  measured  from  the  origin. 
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The  pattern  response  is  given  by 


"a  = Q 


i k X sin  0 


where  Q is  the  source  strength  per  unit  length  and  % is  the  length  of 
the  line.  Integrating,  we  get 


= Q 


i k£/2  sin  6 -i  kJl/2  sin  0 


i k sin  0 


which,  from  an  identity  expressing  the  sine  of  an  angle  in  terms  of  exponentials, 
becomes 


4 kil  , a 

sin  -y  sin  0 

'■a  ■ -U  — 

~Y  sin  0 


, k£ 

1 = -y  sii 


sin  u' 


which  is  just  the  same  as  the  diffraction  of  a slit  in  optics. 

An  important  relationship  for  a plane  rectangular  radiating  surface  is 
obtained  in  the  case  where  the  source  strength  F(x,y)  = f(x)  g(y)  . In 
this  case 


1 k(x  cos  a + y cos  0) 

F(x,y)  e dx  dy 


I f(x) 


i k X cos  a 


iky  cos  6 


dx  g(y)  e 


where  a , 6 , V are  angles  between  direction  of  observation  and  the 


X,  y and  z axes,  respectively. 
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This  very  important  theorem  states  that  if  the  source  strength  for  a 
plane  rectangular  surface  may  be  expressed  as  the  product  of  a function 
f(x)  and  a function  f(y)  , then  the  directional  pattern  is  the  product 
of  two  patterns  in  planes  containing  the  normal  to  the  surface  and  parallel 
to  the  X and  y axes,  respectively.  The  full  significance  of  this  theorem 
will  appear  later  when  we  discuss  rectangular  arrays  of  elements. 

Since  we  have  already  shown  that  patterns  in  the  principal  planes  of 
a rectangular  array  of  point  sources  can  each  be  given  by  a line  of  point 

sources,  we  shall  consider  an  equally  spaced,  in-phase  line  of  sources.  ^ 

The  symbols  that  shall  be  used  for  such  arrays  are  shown  in  Figure  3.  The 
incident  sound  wave  meets  the  line  array  at  the  angle  6 as  shown.  All  | 

phases  shall  be  referenced  to  that  of  the  signal  received  by  an  element  at  \ 

I 

,) 

the  origin  whether  one  is  actually  there  or  not.  Note  the  right-left  \ 

I 

symmetry  of  the  signal  amplitudes. 

Figure  4 is  a representation  of  the  signals  recieved  by  each  of  the 

TT  d 

elements  in  a six  element  line  array  where  u = sin  0 . The  vectors 

above  the  axis  represent  tlie  signals  received  from  the  elements  on  the 
right  and  those  below  the  axis,  the  signals  from  the  left. 

It  is  easily  seen  from  Figure  4 that  the  sum  pattern  of  the  six  element 
line  source  is  given  by,  ^ 

P , = A,  cos  u + A„  cos  3u  + A_  cos  5u  (3)  I 

a6  1 2 3 

i 

and  that  over  the  major  lobe  the  phase  of  the  entire  line  array  is  the 
same  as^that  of  the  reference  signal.  As  0 and  consequently  u increases, 

P^g  will  at  one  point  become  zero  (a  null  In  the  pattern)  and  then  the 
phase  suddenly  reverses  to  180°.  Further  Increases  of  0 and  u produce 


more  nulls  and  phase  reversals. 


The  characteristics  of  P , ana  P.  . are  illustrated  in  Figure  5.  The 

0,0  6,6 

combined  use  of  these  patterns  give  the  necessary  information  for  locating 
underwater  targets.  Note  that  over  a range  of  incident,  angles  near  0°, 
the  difference  between  the  magnitudes  of  the  sum  and  difference  patterns 
increases  nearly  linearly  with  0 . If  the  difference  pattern  leads  the 
sum  pattern  by  90°,  the  target  is  on  the  right  and  if  it  lags  by  90°,  the 
target  is  on  the  left.  In  practice,  the  difference  pattern  is  usually  phase 
shifted  90°  so  that  it  is  in  phase  for  targets  on  one  side  and  out  of  phase 
for  the  other. 

It  must  be  emphasized  that  those  phase  relationships  hold  only  over  a 
relatively  small  range  of  Incident  angles  on  both  sides  of  the  normal.  In 
the  minc^  lobe  regions  we  noted  earlier  that  phase  reversals  occur  as  0 is 
varied  and  where  they  occur  are  not  the  same  for  both  patterns.  Consequently, 
minor  lobes,  unless  drastically  reduced,  will  give  false  information.  For 


-6- 


July  7,  1977 
PMKrcdn 


these  reasons  and  others,  methods  for  reducing  minor  lobes  have  been  developed 
which  are  known  as  tapering  or  shading  the  array.  It  consists  of  a systematic 
variation  of  the  signal  amplitude  coefficients  of  the  elements.  The  process 
will  be  illustrated  by  an  application  to  the  six  element  array  illustrated 
above . 

Note  that  the  pattern  function  was  given  by 

P , = A,  cos  u + A„  cos  3u  + A.  cos  5u  (3) 

uo  1 / J 

which  is  a Fourier  series.  This  suggested  the  use  of  a Fourier  transform 

relating  the  source  strengths  of  the  array  to  the  pattern 

function.  In  one  application,  an  amplitude  distribution  was  derived  to 
provide  the  Gaussian  pattern.  In  this  case,  the  Fourier  transform  is  itself 
a Gaussian  function. 

However,  the  most  common  method  of  tapering  is  that  developed  by 
Dolph  (Reference  2)  and  generally  known  as  the  Dolph-Chebyshev  method  because 
it  sets  the  pattern  function  above  equal  to  an  appropriate  Chebyshev  poly- 
nomial . 

In  Equation  (3)  above,  we  let  x = cos  u and  then  express  the  cosine 
functions  in  powers  of  x using  Chebyshev  polynomials.  Thus,  Equiition  (3) 
becomes 

P , = A,x  + A„(Ax^  - 3x)  + A.(16x^  - 20x^  + 5x) 

O , D 1 Z J 

= 16A^X^  + (4A2  - 20A^)x^  + (A^  - 3A^  + 5A^)x  . (4) 

Tlio  use  of  C'nebyshev  po  1 ynom  I .i  1 s abo’.’e  In  somewhat  trivial  compared  to  the 


really  significant  use  wlilch  is  to  set  the  array  polynomial  above,  equal  to 
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the  Chebyshev  polynomial 


T_(x  x)  = 16x  ^x^  - 20x  ^x^  + 5x  X 

5 o o o o 


(5) 


where  x >1  and  is  given  by  x = cosh{—  cosh  ^r} 


(6) 


where  r = lobe  ratio, 

minor 


The  amplitudes  A^,  and  A^  are  nov;  easily  calculated  by  equating 

coefficients  of  like  powers  of  Equations  (4)  and  (5). 

Figure  6 Illustrates  the  properties  of  the  fifth  degree  Chebyshev 
polynomial  which  applies  to  the  six  element  array.  Over  the  range  of 
-1  ^ x^x  _<  1 , the  Chebyshev  polynomial  function  ranges  between  +1  because 
in  this  range  the  function  is  limited  to  the  cosine  function.  However,  by 
extending  the  independent  variable  x^x  beyond  +1,  the  Chebyshev  polynomials 
can  take  values  to  + infinity.  Thus,  we  can  make  x^x  and  consequently 
T^  (x^x)  as  large  as  possible.  T^  (x^x)  will  then  correspond  to  the  peak  or 
maximum  of  our  response.  On  the  other  hand,  the  minor  lobes  will  all  be 
the  same  and  equal  to  unity. 

This  design  method  is  optimum  in  the  sense  that,  for  a given  minor  lobe 
reduction,  one  obtains  the  narrowest  possible  beamwidths  for  the  equally 
spaced,  in-phase  elements  at  a given  frequency.  We  design  so  that  at 
0 = 90°,  the  Chebyshev  polynomial  is  of  unity  magnitude  just  beyond  the 
last  null. 

There  are  relatively  narrow  limits  on  the  range  of  element  spacing. 


If  the  spacing  is  too  large,  e.g.,  greater  than  X , a second  major  lobe 
would  appear  before  0 = 90°.  Also,  the  method  does  not  work  if  the  spacing 
is  much  less  than  X/2  . For  smaller  spacings,  a slightly  different  meth.od 
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was  Introduced  by  Riblet  (Reference  3)  and  extended  by  Pritchard  (Reference  4). 
In  this  method  which  can  be  applied  to  an  odd  number  of  elements,  a Chebyshev 
polynomial  in  terms  of  an  angle  ip  ~ sin  6 is  used  instead  of  u = ^ sin  0 

This  analysis  leads  to  narrower  patterns  than  those  given  by  the  Dolph-Chebyshev 
technique  and  are  called  super-directive  arrays.  However,  as  the  spacing 
becomes  smaller,  some  of  the  excitation  amplitudes  become  negative,  that  is, 
out  of  phase  with  the  others.  For  very  small  spacing,  the  sum  is  considerably 
less  than  the  sum  of  absolute  values.  Smaller  tolerances  on  the  amplitudes 
are  then  generally  required  and  over-all  efficiency  is  decreased. 

Some  characteristics  of  the  Dolph-Chebyshev  patterns  can  be  determined 
directly  from  the  values  of  x^  , d and  \ . Of  course,  the  latter  depends 
on  the  frequency.  These  characteristics  are  (1)  maximum  allowable  spacing, 
and  (2)  beam  width  at  any  specified  levels  below  the  peak.  These  results 
are  shown  in  Figure  7.  Our  experience  with  transducer  elements  cemented  to 
a rubber  pad  and  mounted  in  a housing  showed  that  since  there  is  a baffle 
effect  that  significantly  reduces  levels  at  incident  angles  beyond  70°,  it 
was  possible  to  design  for  the  last  minor  lobe  to  occur  around  70°.  The 
dashed  line  indicates  the  extension  of  this  tolerance. 

We  shall  now  Illustrate  the  calculations  of  the  excitation  amplitudes 
for  our  six  element  array.  If  we  design  for  a 30  dB  minor  lobe  reduction, 

20  log  r will  be  set  equal  to  30,  which  gives  r = 31.6  . From  Equation  (6) 
we  get  a value  of  1.35  for  x^  . From  Equations  (4)  and  (5)  we  obtain  the 
following  values  for  A . 


^1  = 

15.60 

^2  = 

10.70 

= 4.66 

normalised 

N = 

1 . 000 

'^2  = 

0.6S5 

A. 

J 

= 0.298 
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Figure  8 illustrates  how  these  values  are  used  to  extend  the  line  array 
design  to  that  for  a two  dimensional  array  using  Equation  (2).  Due  to  the 
symmetry,  only  one  quadrant  of  the  array  is  necessary.  Figure  9 gives  the 
normalized  values  for  these  elements.  Frequently  the  outside  corner  element 
is  deleted  so  as  to  confirm  to  a circular  housing.  The  pattern  deterioration 
is  generally  small  but  in  any  case  its  effect  is  easily  calculated.  The 
numbers  in  parentheses  give  one  possible  set  of  turns  on  magnetostrictive 
stacks  or  on  the  transformer  secondaries  used  with  electric  coupled  transducer 
elements. 

The  coefficients  that  provide  an  optimum  sum  pattern  are  unsuitable  for 
an  optimum  difference  pattern.  In  fact,  they  degrade  the  difference  pattern. 
Consequently,  the  difference  pattern  must  be  Independently  optimized. 

With  a slight  modification,  the  Dolph  technique  may  be  employed  to 
optimize  a difference  pattern,  i.e.,  where  the  two  halves  are  in  phase 
opposition.  The  difference  pattern  function, 

Pf  , = A,  sin  u + A-  sin  3u  + A.  sin  5u 

0,6  12  3 

may  be  expanded  in  terms  of  sin  u as  follows: 

Pr  , = A,  sin  u + A„(3  sin  u - 4 sin^  u)  + A. (5  sin  u - 20  sln^  u + 16  sin^ 

0,6  12  3 


Now,  if  we  divide  by  sin  u , we  get. 


-y/- ■ = A,  + A-(3  - 4 sln^  u)  + A_(5  - 20  sin^  u + 16  sin^  u) 

sin  u 1 2 3 


2 2 

If  wo  substitute  1 - cos  u for  sin  u and  x for  cos  u , the  expression 
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above  becomes , 

-7-^  = (A,  - A-  + A.)  + (4A-  - 12A  ) + 16A,x^  . 

sin  u 1 Z i Li  i 

We  now  equate  this  relationship  to  a 4th  degree  Chebyshev  polynomial  and 

proceed  as  before  with  the  sum  pattern.  Please  recall,  however,  that  we 

used  a 5th  degree  polynomial  for  the  sum  pattern.  Also,  note  that  since 

P.  , = sin  u T, (x  x)  , it  Is  zero  at  0=0  and  then  the  absolute  value 
0,6  4 o 

Increases  rapidly  for  moderately  small  values  of  9 on  both  sides  of  zero. 

The  use  of  Equation  (2)  to  extend  the  line  array  to  a plane  is  shown 
In  Figure  10.  Here,  the  B factors  are  just  the  excitation  coefficients 
for  an  optimum  sum  pattern.  In  fact,  if  all  elements  were  connected  in 
phase,  we  would  obtain  an  optimum  vertical  sum  pattern,  but  when  the  right 
and  left  halves  are  connected  in  phase  opposition,  we  obtain  the  optimum 
horizontal  difference  pattern. 

Figure  11  gives  relative  source  amplitudes  for  the  optimum  difference 
pattern.  The  distinct  difference  between  these  coefficients  and  those  for 
the  sum  pattern  is  quite  evident  by  noting  that  the  elements  in  the  second 
column  now  have  the  greatest  excitation  instead  of  the  first  column  as  in 
the  case  of  the  sum  pattern. 

No  direct  method  of  finding  x^  was  available  for  the  difference  pattern 
but  minor  lobe  reductions  could  be  calculated  for  assumed  values  of  x 

o 

(see  Reference  1)  and  plotted  as  shown  in  Figure  12.  Thus,  these  plots 
permit  us  to  select  the  appropriate  x^  for  a given  minor  lobe  reduction. 

Figure  5 presents  a linear  plot  of  the  optimum  sum  and  difference 
patterns.  When  this  method  for  the  optimum  difference  pattern  was  developed, 
it  appeared  that  it  would  apply  only  to  an  even  number  of  array  columns. 
However,  Geoffrey  Wilson  (Reference  5)  has  shown  that  it  can  bo  applied  to 
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an  odd  number  of  columns,  in  which  case  the  central  column  simply  is  not 
used  to  form  the  difference  pattern. 

Figure  13  is  of  academic  interest.  It  represents  the  locus  of  the 
resultant  signal  from  1/2  of  a transducer  designed  for  (a)  an  optimum  sum 
pattern  and  (b)  an  optimum  difference  pattern.  It  is  the  resultant  of  the 
vectors  shown  in  Figure  4.  The  sum  pattern  response  is  the  projection  on 
the  x-axis  and  that  of  the  difference  pattern  is  the  projection  on  the  y-axis. 
It  clearly  shows  why  the  minor  lobes  are  so  small. 

Tilted  beams  and  end-fire  arrays  are  obtained  by  introducing  delay  lines 
or  phase  shift  networks  into  the  circuits  of  each  element  or  maybe  column 
or  row  of  the  array  so  that  for  a given  specified  direction,  the  outputs 
of  signals  (including  delays  or  phase  shifts)  will  be  in  phase.  The 
tapering  procedures  are  generally  similar  to  those  described  earlier. 

Most  of  the  development  for  end-fire  arrays  (beam  tilted  to  90°) 
resulted  in  super directive  arrays  with  some  characteristics  similar  to 
those  of  Riblet  (Reference  3)  and  Pritchard  (Reference  4). 

A somewhat  different  approach  to  the  production  of  tilted  sum  and 
difference  patterns  was  developed  by  W.  J.  Hughes  and  W.  Thompson,  Jr. 
(Reference  6).  1 shall  demonstrate  its  use  by  an  application  to  a tilted 
sum  pattern  which  may  be  written, 

P^  = Aj  cos(u  - (|))  + A2  cos  3(u  - ()))  + cos  5(u  - (J))  + . . . 

where  ip  is  the  required  phase  delay  for  elements  on  one  side  of  center 
and  ph*e  advance  for  elements  on  the  other  side.  Using  the  trigonometric 
identity  for  the  cosine  of  the  difference  of  two  angles,  the  equation  above 
becomes , 
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= [Aj^  cos  4>  cos  u + cos  3(j)  cos  3u  + A^  cos  5<p  cos  5u  + . . • ] 
+ [A  sin  (p  sin  u + A^  sin  34>  sin  3u  + A^  sin  5(J)  sin  5u  + . . . ] 

Now  for  a given  tilt  angle  <p  , the  sines  and  cosines  involving  <{)  are 
constants.  Therefore,  the  pattern  function  becomes 


cos 

U + A^ 

COS 

3u 

+ 

^3 

cos 

5u 

+ . 

. .] 

+ [a;; 

sin 

u + A^’ 

sin 

3u 

+ 

"3 

sin 

5u 

+ . 

. .] 

where  A*  = A cos  i , and  A"  = A sin  d)  , etc. 
n n n n 

We  note  that  P is  now  the  sum  of  a sum  pattern  and  a difference 
a ^ 

pattern.  Since  the  latter  is  in  phase  quadrature  with  the  former,  it  must 

be  shifted  90°  before  combining  the  two.  Although  the  method  is  applicable 

to  any  values  of  A , tapering  is  easily  achieved.  The  unpriraed  A 
m n 

coefficients  are  determined  as  before  by  using  Chebyshev  polynomials.  The 
Hughes  - Thompson  technique  is  also  applicable  to  tilted,  difference  patterns 
and  end  fire  arrays.  It  is  Interesting  to  note  that  all  pattern  tapering 
techniques  discussed  so  far  have  used  Chebyshev  polynomials. 

An  array  tapering  for  an  entirely  different  purpose  was  developed  by 
W.  J.  Trott  (Reference  7).  Its  purpose  was  to  provide  plane  waves  in  a 
region  of  the  near  field  of  a relatively  large  array  in  order  to  calibrate 
transducers  in  the  near  field. 

It  is  well  known  that  the  Intensity  along  the  axis  of  a circular 
radiating  piston  whose  diameter  is  large  compared  to  the  wavelength  has 
numcroi^  maxima  and  minima  (actually  nulls)  in  the  near  field.  It  is  also 
known  that  the  total  number  of  side  lobes  appearing  in  the  far  field  is 
just  twice  the  number  of  maxima  which  occur  directly  in  front  of  the  piston. 
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The  above  remarks  suggest  that  if  minor  lobes  are  eliminated  or  greatly  reduced, 
this  axial  variation  should  likewise  be  eliminated  or  reduced.  Indeed  this 
is  the  case.  Applying  these  principles,  Trott  adjusted  the  excitation 
amplitudes  of  a large  array  of  elements  in  order  to  provide  an  essentially 
plane  wave  region  directly  in  front  of  the  array.  In  his  example,  a 100  cm 
diameter  array  provided  a cylindrical  space  in  front  of  the  array  about  50  cm 
in  diameter  and  about  the  same  in  length  wherein  there  were  plane  waves. 

Although  it  was  not  expressly  stated,  all  the  previous  discussion  was 
based  on  linear  acoustic  waves.  When  very  intense,  so-called  finite  sound 
waves  are  projected,  non-linear  or  distorted  waves  are  produced  with  a whole 
host  of  interesting  properties.  For  example,  if  two  pure  tones  are  projected 
into  the  same  medium,  harmonics  of  the  waves  are  produced,  as  well  as  sum 
and  difference  frequency  waves.  The  highly  directive  character  of  these 
difference  frequency  beams  was  first  recognized  by  Peter  J.  Westervelt 
(Reference  8 ) . 

In  the  following  discussion  we  are  primarily  interested  in  the  difference 
frequency  that  develops  in  the  projected  acoustic  field.  Usually  the 
frequencies  of  the  high  intensity  primary  waves  are  close  together  so  that 
the  difference  frequency  is  much  smaller,  say  one-tenth  or  less. 

When  the  near  field  primary  wave  absorption  is  quite  small,  the  inter- 
action zone  will  extend  beyond  the  mean  Rayleigh  distance  r^  . The  inter- 
action zone  is  that  region  in  front  of  the  transducer  face  where  wave 
distortion  occurs  and  acoustical  energy  is  pumped  from  the  primary  waves 
into  the  distortion  waves,  including  the  difference  frequency  wave.  Since 
this  latter  is  of  much  lower  frequency  than  any  of  the  other  waves  and 
consequently  attenuated  rmuh  less  than  the  others,  it  usually  will  be 


transmitted  to  greater  r.inges  tli.ut  the  others  even  though  its  equivalenr 
field  pressure  Is  much  below  that  of  the  primary  waves.  The  so-called 
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pumping  occurs  until  the  spreading  loss  of  the  primary  waves  reduce  the 
Intensity  to  such  a level  that  distortion  no  longer  exists.  Since  the 
virtual  sources  that  produce  the  difference  frequencies  are  essentially 
limited  to  the  main  beams  of  the  primary  waves,  the  low  frequency  difference 
pattern  will  be  essentially  the  same  as  those  of  the  much  higher  frequency  primary 
patterns  or  even  a little  sharper.  In  fact,  its  far  field  directivity 
response  is  essentially  equal  to  the  product  of  the  directivity  responses 
of  the  primary  patterns.  One  very  Important  consequence  of  this  is  the 
almost  total  absence  of  minor  lobes  in  the  difference  frequency  pattern. 

In  order  to  achieve  the  same  directivity  with  a conventional  source 
operating  at  the  difference  frequency,  the  source  diameter  would  have  to 

exceed  the  diameter  of  the  parametric  array  by  the  ratio  of  the  mean  | 

I : 

primary  frequency  to  the  difference  frequency. 

The  bandwidth  of  the  difference  frequencies  may  be  very  broad  because  i 

it  is  really  determined  by  the  Q of  the  transducer  that  projects  the 
primary  waves. 

Somewhat  similarly,  difference  frequency  waves  can  be  produced  which 
have  a constant  beam  width  over  very  wide  frequency  bands  (e.g.,  greater 
than  two  octaves).  Of  course,  sum,  difference  and  tilted  beams  can  be 
produced . 

When  the  absorption  coefficient  of  the  primary  waves  is  large,  the 
attenuation  may  be  so  great  that  all  of  the  Interaction  may  take  place 
within  the  near-field  of  the  transducer,  l.e.,  at  ranges  less  than  r 

o 

In  this  case,  the  interaction  is  essentially  confined  to  a cylinder  whose 
dlamet£W^  is  that  of  the  transducer  radiating  face  and  of  length  roughly 
determined  by  the  absorption  coefficient.  This  interaction  zone  is  really 
an  extension  of  the  source  itself.  Indeed,  the  beam-width  of  this  virtu.-il 


end-fire  array  is  given  by  the  same  relationship  as  that  for  an  end-fire 
array  and  is  always  less  than  that  for  the  "spreadlng-loss-limited"  case. 

Parametric  receiving  arrays  are  formed  in  a fluid  by  projecting  a 
finite-amplitude  "pump  wave"  into  the  medium  to  serve  as  a "carrier"  for 
the  weak  incoming  signal  whose  frequency  equals  the  difference  frequency 
of  the  transmitted  waves. 
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